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Black hole quasinormal frequencies are complex numbers that encode information on how a black 
hole relaxes after it has been perturbed and depend on the features of the geometry and on the 
type of perturbations. On the one hand, the examples studied so far in the literature focused on the 
case of black hole geometries with singularities in their interior. On the other hand, it is expected 
that quantum or classical modifications of general relativity may correct the pathological singular 
behavior of classical black hole solutions. Despite the fact that we do not have at hand a complete 
theory of quantum gravity, regular black hole solutions can be constructed by coupling gravity to an 
external form of matter, sometimes modeled by one form or another of nonlinear electrodynamics. 
It is therefore relevant to compute quasinormal frequencies for these regular solutions and see how 
differently, from the ordinary ones, regular black holes ring. In this paper, we take a step in this 
direction and, by computing the quasinormal frequencies, study the quasinormal modes of neutral 
and charged scalar field perturbations on regular black hole backgrounds in a variety of models. 



I. INTRODUCTION 

The problem of understanding how avoidance of singu- 
larities may be possible in general relativity goes beyond 
formal importance and is not at all new. In 1968 Bardeen 
constructed the first example of a regular black hole, i.e., 
a regular, nonsingular geometry with an event horizon, 
satisfying the weak energy conditions [Tj . Despite its the- 
oretical relevance, Bardeen's regular black hole solution 
lacked, for several years, a satisfactory physical interpre- 
tation. The reason is that it is not a vacuum solution 
of Einstein's equations and, in order to generate it, it is 
necessary to introduce some external form of matter or 
a modification to gravity. The solution was obtained by 
introducing an ad hoc energy momentum tensor, regu- 
lar and bounded, decaying at infinity, and satisfying the 
weak energy conditions. However, no fundamental physi- 
cal motivation for this choice had been given, until Ayon- 
Beato and Garcia reobtained the solution by describing 
it as the gravitational field of some sort of nonlinear mag- 
netic monopole [2]. 

Bardeen's work has been followed by several other ex- 
amples, motivating deeper analyses of how singularity 
avoidance may be possible in general. Other solutions 
have been proposed in the literature and some attention 
has been directed to theories of gravity coupled to non- 
linear electrodynamics. Solutions have been discussed in 
different contexts, and a few of interest to us are, in ad- 
dition to Refs. [TJ those analyzed in Refs. [3HH] ■ It 
is worth mentioning that an important viable example 
of a black hole with a regular center was constructed by 
Dymnikova, with a de Sitter core smoothly connecting 
to a Schwarzschild outer geometry [9]. Analysis of reg- 
ular black hole solutions continued in several directions, 
see, e.g., Refs. [TU] [TT] . Another important step toward 



understanding the absence of singularities in general rel- 
ativity was taken by Borde who showed that, for a large 
class of black hole solutions, absence of singularities was 
related to a change in the topology beyond the event hori- 
zon [T2] . Borde's arguments clearly demonstrated the im- 
possibility of proving general singularity theorems when 
the strong energy condition or the existence of global hy- 
perbolicity were not assumed. 

Now, the interior is a hidden region by definition. The 
region that connects the interior to the exterior is the 
horizon and its near horizon region. Thus, a way to peer 
into the interior is to perturb the horizon. More im- 
portant, quantum processes occur in its neighborhood, 
giving a glimpse of phenomena that unites quantum me- 
chanics and gravitation. Indeed, the Hawking radiation 
has its origin in the vicinity of the horizon and the black 
hole entropy is believed to come out from degrees of free- 
dom located at the horizon. Thus, the horizon is the 
region that may give a glimpse of what the black hole 
interior is. By poking the horizon something from the 
inside might pop up. 

One way to poke the horizon is to perturb space- 
time, create quasinormal modes (QNMs), and analyze 
the result of the perturbation. Quasinormal frequencies 
(QNFs) are complex numbers that encode information 
on the system's relevant parameters and on its relaxation 
after it has been perturbed. Quasinormal modes are re- 
lated to the classical evolution of a system. They can 
also reveal instability. A way to study QNMs is through 
a WKB approximation as was done first in [15] . This 
was further developed by Iyer and Will [H] [TS] to study 
the Schwarzschild black hole and by [TBT - ET] to study 
perturbation in the Reissner-Nordstrom case as well as 
charged scalar perturbations. Modes with large imag- 
inary parts were analyzed, in connection with a possi- 
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ble horizon area quantization law, in Schwarzschild and 
Reissner-Nordstrom see [22j [23] . For reviews see [24U26] . 

All this work has been done for vacuum black holes. 
It is important to proceed for regular black holes. A 
first step was made by Fernando and Correa [27] who 
computed the QNFs for Bardeen's solution [HE]. I n 
addition in Ref. 28 a study of the QNMs for the solution 
presented in [8] was performed. In this paper, we intend 
to study QNFs for the solutions given in Refs. [THH1 - 

II. SOME REGULAR BLACK HOLES 

In the above context of gravity coupled to some form 
of matter, regular solutions were obtained from a proto- 
typical action of the form 



(1) 



where g is the determinant of the metric g^, R is the 
scalar curvature, and Jz? represents the Lagrangian of 
the matter fields. For the case of nonlinear electrody- 
namics then J>f — Jz? (F) is a nonlinear function of the 
electromagnetic field strength with F — 
The general line element 
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when presented for spherically symmetric regular black 
hole solutions takes the form 



ds 2 



dv 

f {r)dt 2 + — +r 2 (d9 2 + sm 2 9 d(j> 2 ) , (3) 
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where (t, r, <fi) are the usual space-time spherical coor- 
dinates and the lapse function / = f(r) depends on the 
specific form of underlying matter. 

In [T] the function takes a particularly simple form 
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with a = const and m being the mass of the solution. 
This implies a specific matter energy-momentum tensor 
that is de Sitter at the core and vanishes away at infinity 
as a magnetically charged solution in the context of a 
specific nonlinear electrodynamics with Lagrangian given 
by JS? = (3/(2sa 2 ))( %/2a 2 F/(l + %/2a 2 F)) 5 / 2 , where s = 
\a\/2m and a is the magnetic charge [2]. Depending on 
the relative values of m and a, Eq. (§ can have two, one 
or zero horizons. 

In [3] the function also takes a simple form 
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which is a variant from the original Bardeen's proposal 
PQ with a = constant. This also implies a specific matter 
energy-momentum tensor that is de Sitter at the core 
and vanishes away at infinity. As well, depending on the 
relative values of m and a equation ^ can have two, one 
or zero horizons [3]. 



The basic properties of the Lagrangian _£f (F) leading 
to electrically charged configurations have been discussed 
originally in Bronnikov 4, 5 and used by Dymnikova [6 
to show that regular electrically charged solutions com- 
patible with the weak energy condition require a de Sitter 
center and to construct a prototypical example. Specifi- 
cally, Dymnikova's solution is obtained from a nonlinear 
electrodynamic theory with a Hamiltonian-like function 

(see [6] for details) of the form Jrf? — P (l + ay—P) , 
cwhere = _Sf (F) F^ {P = P^P^) leading to a 
solution of the above form with 
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The parameter ro in the solution Q is a length scale re- 
lated to the total mass m and the charge q by the relation 
ro = irq 2 /(8m). The above solution reproduces asymp- 
totically, for r — > oo, the Reissner-Nordstrom behavior, 
while near the center, for r <C the solution approx- 
imates to de Sitter. Depending on the value of r , the 
above regular solution may present two distinct, one, or 
no horizons (see Fig. [T]) . Another electrically charged 
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Figure 1. The behavior of the metric function f(r) given in 
Q for different values of the parameter ro. The thick dashed 
line refers to the case where the two horizons overlap. The 
bottom-right inset represents the values of the outer (inner) 
horizon r+ (r_) vs the parameter ro. The mass m has been 
normalized to unity and ro varied up to its critical value above 
which the horizons disappear and the geometry is regular. 

solution discussed in the literature is the one presented 
by Ayon-Beato and Garcia in Ref. [7J. It takes the form 
([3]) with / given by 



/ = 1 
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and it is obtained from a nonlinear electrodynam- 
ics with Lagrangian density _S? = df^^fz^y^ 



. where X = y— 2q 2 F, and m and q are 



3m X 5 ^ 2 (3-2X) 
2g3 (1 _x) 7 /2 

associated with mass and charge, respectively. The above 
solution is asymptotically flat and behaves like Reissner- 
Nordstrom, at leading order, asymptotically. Depending 
on the values of the charge and mass, as Dymnikova's so- 
lution ([6]), it has two distinct inner and event horizons, for 
values of the charge smaller than a critical value q c , two 
degenerate horizons for q = q c , and becomes a globally 
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regular geometry for q > q c . This solution, as discussed 
in Ref . [4] , presents a problem related to the presence of 
cusps in the electromagnetic Lagrangian. As pointed out 
in Ref. [4 , such a solution should be taken with care, and 
here we will use it as a useful working example and to 
check our numerics against known results. 

Amongst the various models, the one of Ref. [5] is 
slightly more elaborate and the solution arises from the 
system of gravity coupled to a phantom scalar field. 

All the solutions considered in this paper are summa- 
rized in Table I. 
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Table I. The models of Refs. [18] were constructed in such 
a way to have at least Reissner-Nordstrom asymptotics (i.e., 
order of 1 + 1/r + l/r a with a > 2, see Ref. [1,7]) or at 
least Reissner-Nordstrom-de Sitter asymptotics (i.e., order of 
r 2 + l/r + l/r a with a > 2, see Ref. [8] for the case a = 3). 
The solution [TJ [2] can be put in a nonlinear electrodynam- 
ics framework. The solution given in [3] is a minimal regular 
solution with the parameter m measuring the mass and the 
parameter a measuring the deviation of the nonsingular so- 
lution from Schwarzschild, which is reproduced for a = 0. 
The models of Refs. 4- 7 were constructed in the context of 
nonlinear electrodynamics coupled to gravity using different 
functional forms for . The parameter ro in [2H6] is a length 
scale related to the electric charge and q in [7] is the electric 
charge itself. The solution proposed in [8] was constructed 
coupling gravity to a phantom scalar field with a potential 
U(b, c, po) where the parameters b, c, and po characterize the 
features of the potential. All the details can be found in the 
original references. 



III. QUASINORMAL FREQUENCIES 

The computation of QNMs and their QNFs has been at 
the center of great attention both for the astrophysical 
relevance in relation to gravitational wave observation 
and for its formal importance. Extensive computations of 
QNFs for various black hole geometries are now available 
and most results have been reviewed [TBH^tJ] . 

Here we extend the computation of the QNFs to regu- 
lar black hole geometries; see also [27] [35] . Although we 
will consider the specific solutions listed in Table I, it is 
worth noticing that the QNFs for uncharged scalar per- 
turbations, as well as for tensor perturbations, depend on 
the specific form of the matter Lagrangian used to obtain 
the background solution only through the explicit func- 
tional form of the lapse function /. The same is also true 
for the case of charged scalar and vector perturbations, 



as long as the weak held limit of the nonlinear electrody- 
namics reproduces the standard Maxwell's theory. 

A. Neutral scalar perturbations 

The computation of the QNFs for regular black hole 
geometries proceeds as in the well known cases, and this 
section will be devoted to those arising from neutral 
scalar field perturbations. The equation for these per- 
turbations takes the usual form, 
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d„ (V99^d v ) = 
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where g is the determinant of the metric tensor given 
in ([3]) and the lapse function / depends on the model con- 
sidered. Equation ([8]) can be separated by decomposing 
the scalar perturbation into appropriate harmonics, 



It, 



(9) 



and by introducing the tortoise coordinate, 

dx = dr/ f . 

One then can rewrite the radial part of 
Schrodinger form, 
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with V given by 



v = f 



\ r 
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and where the indices I, m have been suppressed from 
ipim,. The boundary conditions, appropriate for the com- 
putation of the QNFs, must force the solution near the 
event horizon, x = —00 (at infinity, x — +00), not to 
generate outgoing (ingoing) waves. These can be written 
as 

ip ~ e ~ ZUJX , x -> -00 , 
<p ro e +iux , x -> +00 . 

In the present case, most available methods to compute 
QNFs would work. 

Here, we will adopt the most direct way that makes use 
of WKB approximation, originally discussed in Ref. [2] , 
and that does not require any special modification. Re- 
sults for the quasinormal frequencies for neutral scalar 
perturbations are tabulated, in Appendix A, in Table II 
for the model given in Refs. [U [2], in Table III for the 
model given in Ref. i3] , in Table IV for the model given in 
Refs. [HEO, hi Table V for the model given in Ref. (B], in 
Table VI for the model given in Ref. [7] , and in Table VII 
for the model given in Ref. [8] . The tables list the values 
of the frequencies at third order in the WKB approxima- 
tion, as results to this order allow a direct comparison 
with the results presented in [Kl HH] ■ One can go up to 
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sixth order in the WKB approximation and this allows 
to check the convergence of the approximation and how 
the accuracy improves going to higher order. 

Figure 2 illustrates the real and imaginary parts of 
the frequencies for the solution of Dymnikova [S] and for 
sample values of the parameters, multipole number I, and 
overtone number n, up to sixth order WKB, thus allowing 
to test the convergence of the approximation. 
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Figure 2. The real (upper panel) and imaginary (lower panel) 
parts of the QNFs from neutral scalar perturbations for sam- 
ple values of the overtone and multipole numbers for the 
model of Dymnikova 6 when higher order terms in the WKB 
approximation are included in the computation. The symbols 
o (+) refer to ro = 0.2 (ro = 0.3) and we have normalized m 
to unity. The overtone and multipole numbers are indicated 
in the figure. 

Because of its connection with a possible black hole 
area quantization law, it seems interesting to study the 
behavior of the QNFs in the limit of the large imagi- 
nary part. The problem is rather similar to the case 
of charged black hole solutions studied, for example, in 
Ref. |23| . This comes as no surprise, since the regular 
black hole solutions considered here have, at leading or- 
der, an asymptotic structure at least of the Reissner- 
Nordstrom type or Reissner-Nordstrom-de Sitter type, 
i.e., order of 1 + 1/r + l/r a with a > 2 (see [1,7]), or 
order of r 2 + 1/r + l/r a with a > 2 (see [8] for the 
case a — 3), respectively. Although formally the process 
of monodromy matching can be performed in the same 
way, one can see that in the limit r — > 0, the leading term 



in the potential (11) behaves as 1(1 + l)r~ 2 , rather than 
r~ 6 . This suggests that the QNFs, in the limit of the 
large imaginary part will depend on the multipole num- 
ber Following the same steps as in [23], one finds that 
asymptotically the frequencies will satisfy the following 
relation: 



e>~ = -(1 + 2 cos(7iVl + J(Z + l))) 

-e^(2 + 2 cos(vrv/l + 1(1 + 1))) , (12) 

where (3 and /3_ are the inverse temperature at the outer 
and inner horizons. The above formula suggests that the 
real part of the QNFs does not asymptote to a constant 
value, nor display any periodic behavior. It may be inter- 
esting to study what happens in the extremal case. If the 
analogy with the Reissner-Nordstrom black hole is valid 
also in this case, then, as the charges approach the ex- 
tremal value, the asymptotic value of the real part of the 
QNFs would spiral toward the asymptotic Schwarzschild 
QNFs [Tni [2D] - To test this, however, a more suitable nu- 
merical approach is necessary and this is left for future 
work. 

B. Charged scalar perturbations 

The computation of QNFs can be extended to charged 
scalar perturbations. The problem is similar to the 
Reissner-Nordstrom case treated in Refs. [T6HT9] where 
QNFs associated with charged scalar perturbations were 
computed. In particular, Ref. [T5] followed the analysis 
of Ref. [17] that studied the evolution of scalar pertur- 
bations around collapsing charged black holes, indicating 
that scalar perturbations are radiated away, at a slower 
rate than the neutral ones, and follow an inverse power- 
law behavior at the future outer horizon, while displaying 
a decaying behavior accompanied by oscillations at the 
future outer horizon. 

The general equation for the (complex) scalar pertur- 
bations can be written formally in the same way as [17] . 

^ ab g ab - ieA a g ab (24>, b - ieA b 4>) - ieA a , b g ab (b = 0, (13) 

with e being the (constant) charge of the scalar field. 
Since we are concerned with perturbations, the electro- 
magnetic potential for the black hole is determined by 
the ordinary Maxwell's theory and can be written, up to 
an additive integration constant c, as 
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dt + c 



The integration constant is to be determined by regular- 
ity at the horizon for the gauge field. In fact, as far as the 
evaluation of the QNFs is concerned, its precise value is 
unessential since it would only produce a uniform shift in 
the real part of the frequencies. In the following, we will 
set c = 0. Equation ( 13 ) can be separated by decompos- 



ing the scalar perturbation into appropriate harmonics, 
and using tortoise coordinates it is possible to rewrite the 
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above equation in Schrodinger form, 
d 2 



+ V{x) - {w + eq/r)' 
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(14) 



where V is given by (11). The boundary conditions rel- 



evant for the QNFs computation are 

,-i(u+eq/r+)x 
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The WKB method works in this case too and the com- 
putation traces over the Reissner-Nordstrom case, with 
the difference being the lapse function f(r), requiring a 
small adaptation in the numerics. In the case of charged 
perturbation some more work is necessary, but the same 
approach of analytic continuation used in Ref. |18j can 
be adopted here: first fix all the parameters including 
the multipole and overtone numbers, then maximize the 
potential as a function of the radial distance, and finally 
find the value of omega that satisfies the WKB relation 
for the QNFs as a numerical function analytically contin- 
ued into the complex domain (the potential is generally a 
complex function). The results for the QNFs for charged 
scalar perturbations are listed in Appendix B, in Table 
VIII for the model given in Refs. pQ [2], in Table IX for 
the model in Ref. [3], and in Table X for the model in 
Ref. 0. 

Figures [3] and [3] illustrate the behavior of the QNFs 
with respect to the charge q for the solution of [6j. It is 
possible to notice that for vanishing q the solution repro- 
duces the Schwarschild geometry. In this limit, scalar 
perturbations become uncharged and the value of the 
QNFs must be smoothly connected with the QNFs for 
uncharged scalar perturbations on Schwarzschild black 
holes, indicated by the black dots in Figs. 3]|T 
[5] and [6] illustrate the QNFs for the I 
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mode for different values of the charge of the perturba- 
tion in the Reissner-Nordstrom case and in the regular 
black hole case. The choice of these specific values for 
the multipole and overtone numbers was motivated by 
the desire to compare our results with those of Ref. [TH] • 
One notices (upper panel of Fig. [6]) that the real part of 
the QNFs follows a behavior analogous to the Reissner- 
Nordstrom case; i.e., Rw grows with the charge and it 
is larger for charged perturbations than for neutral ones. 
This is a natural consequence of the fact that the regular 
black hole geometries considered asymptote, at leading 
order, to Reissner-Nordstrom. This indicates that, in the 
regular case as well, the late time behavior of the quasi- 
normal ringing will be dominated by the neutral scalar 
perturbations. The lower panel illustrates the imaginary 
part of the QNFs for Reissner-Nordstrom and for the 
regular black hole solution (J3j|. While for small values 
of the charge q, i.e., in the region where the deviations 
from the two geometries are small, the values of 9 uj are 
basically superposed, the discrepancy increases with the 
charge signifying that in the large q region the exponen- 
tial damping of these modes will occur less prominently 
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Figure 3. The real (o) and imaginary (+) parts of the QNFs 
with I = 1 from charged scalar perturbations for sample val- 
ues of the overtone and multipole numbers (as indicated in 
the panels) for the model of Dymnikova [6|. The dark cir- 
cles represent the values of the QNFs for the Schwarzschild 
case that are continuously connected with the values for the 
Dymnikova solution in the limit of zero charge (q = 0). 
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Figure 4. Same as Fig. [3]with I = 2. 



than in the Reissner-Nordstrom case. The coincidence in 
the imaginary parts of the QNFs is observed in the regu- 
lar case as well, but the feature of the profile, namely the 
fact that 9 u) undergoes a transient increase before drop- 
ping for large values of the charge, is much less prominent 
in the regular case, where the increase in 3w for moder- 
ate values of q is small. 

IV. CONCLUSIONS 

In this work we have computed the neutral and charged 
scalar QNFs for regular black hole geometries. The pro- 
totype example of such regular solutions is the one pro- 
posed by Bardeen many years ago Q] (see also [2]), and 
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being perturbed can shed light on some fundamental as- 
pects of quantum gravity. In turn, knowledge of quantum 
gravity should provide us with a better understanding 
of black holes and eventually suggest a possible resolu- 
tion of the singularity problem. In fact, the presence of 
singularities certainly signals a limitation of our under- 
standing, if not a breakdown, of general relativity. This 
pathological behavior is usually believed to disappear in 
a full theory of quantum gravity that would provide a 
consistent framework to test the well known semiclassi- 
cal arguments predicting the evaporation of black holes. 



Figure 5. The behavior of the real part of the QNFs for / = 3 
and n = in the Reissner-Nordstrom case and in the regular 
black hole case for the model given in 0, see Eq. 0. The 
symbols refer to q = (+), 0.1 (jj), 0.3 (*) for the regular case, 
and q = 0(b), 0.1(1)), 0.3(A) for the Reissner-Nordstrom 
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Figure 6. Same as Fig. [5]for the imaginary part of the QNFs. 



several examples were constructed following this original 
example. Such solutions are spherically symmetric, mul- 
tihorizon, and singularity-free and display an asymptotic 
behavior analogous, at leading order, to at least Reissner- 
Nordstrom or Reissner-Nordstrom-de Sitter, while hav- 
ing a near-horizon behavior similar to Schwarzschild. 
The WKB method can be implemented in the present 
case, and it has been adopted to compute the QNFs for 
the variety of models of Refs. [THE]. The results have 
been tested in the limit of Schwarzschild and Reissner- 
Nordstrom black holes reproducing the known results of 
Refs. [151 HE] when the same order in the WKB approx- 
imation is used and the convergence has been tested up 
to sixth order in the WKB expansion. The QNMs of 
Bardeen's model have also been studied in [UJ HE] • 

Aside of the generalization to other types of perturba- 
tions, the use of more sophisticated numerical methods, 
like those used in Ref. [19] . may be useful to check the 
accuracy of the WKB approximation as well we to test 
the asymptotic behavior of the QNFs. These problems 
are left for future work. 

Certainly, knowledge of how a black hole rings after 



7 



Appendix A: Tables of QNFs for neutral scalar perturbations 
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n i 

U.l 
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n q 
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U.lU4y — 


r, n 1 1 /in 


U.llWU 


« n 1 1 op: 


U.lUyo — 


2 U.lUlo 




I 


n nsQ5 — 


L yj . OtJ^t-L 


0.0921 




0866 — 


j 391 Q 


1 





0.2916- 


i 0.0978 


0.2959 


- % 0.0967 


0.3132- 


i 0.0906 




1 


0.2629- 


i 0.3069 


0.2686 


- i 0.3028 


0.2886 - 


« 0.2820 




2 


0.2245 - 


i 0.5259 


0.2323 


— i 0.5187 


0.2548 - 


i 0.4831 




3 


0.1750- 


i 0.7474 


0.1860 


- i 0.7368 


0.2125 - 


2 0.6869 


2 





0.4840 - 


i 0.0966 


0.4909 


- i 0.0956 


0.5191 - 


« 0.0901 




1 


0.4641 - 


j 0.2954 


0.4721 


- % 0.2920 


0.5034 - 


i 0.2740 




2 


0.4328 - 


i 0.5027 


0.4424 


- i 0.4966 


0.4777- 


i 0.4645 




3 


0.3940 - 


I 0.7148 


0.4058 


- i 0.7059 


0.4454 - 


i 0.6595 



Table II. QNFs for neutral scalar perturbations for q — 0.1, 0.3, 0.6 and m = 1 for the model of Refs. [T] [2]. The model 
reproduces the original example constructed by Bardeen using a nonlinear electromagnetic theory that displays, in the weak 
field limit, a stronger behavior as compared to the ordinary Maxwell's one. In this model the gravitational field can be 
interpreted as that of a nonlinear magnetic monopole. 



1 


n 


a 


= 0.1 


a = 


: 0.4 


a = 


0.7 








0.1046 


- i 0.1148 


0.1034- 


i 0.1093 


0.0923 - 


i 0.0986 




1 


0.0891 


- i 0.3541 


0.0835 - 


i 0.3424 


0.0562 - 


i 0.3281 


1 





0.2913 


- i 0.0978 


0.2946 - 


i 0.0948 


0.3018- 


i 0.0854 




1 


0.2624 


- i 0.3068 


0.2655 - 


i 0.2973 


0.2649 - 


i 0.2691 




2 


0.2238 


- i 0.5258 


0.2264- 


i 0.5104 


0.2108- 


i 0.4677 




3 


0.1741 


- i 0.7473 


0.1764- 


i 0.7262 


0.1418- 


i 0.6727 


2 





0.4835 


- i 0.0966 


0.4892 - 


i 0.0940 


0.5037- 


« 0.0855 




1 


0.4636 


- i 0.2953 


0.4700 - 


i 0.2870 


0.4823 - 


i 0.2601 




2 


0.4321 


- i 0.5025 


0.4392 - 


i 0.4881 


0.4445 - 


z 0.4422 




3 


0.3931 


- % 0.7146 


0.4009 - 


i 0.6942 


0.3946 - 


2 0.6310 



Table III. QNFs for neutral scalar perturbations a = 0.1, 0.4, 0.6 and m — 1 for the solution of Refs. [3]. 
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n 


ro = 


= 0.1 


ro = 


= 0.3 


ro = 


: 0.4 








0.1091 - 


i 0.1151 


0.1189- 


/ 0.1119 


0.1200- 


i 0.1065 




1 


0.0940 - 


% 0.3550 


0.1021 - 


i 0.3482 


0.0940 - 


i 0.3401 


1 





0.3019- 


i 0.0988 


0.3303 - 


% 0.0990 


0.3503- 


2 0.0962 




1 


0.2742 - 


i 0.3094 


0.3058 - 


i 0.3080 


0.3252 - 


2 0.2980 




2 


0.2374 - 


i 0.5297 


0.2729 - 


i 0.5264 


0.2897- 


i 0.5098 




3 


0.1902 - 


% 0.7523 


0.2319 - 


i 0.7472 


0.2451 - 


i 0.7249 


2 





0.5007- 


i 0.0977 


0.5472 - 


% 0.0982 


0.5810- 


2 0.0957 




1 


0.4815 - 


i 0.2984 


0.5306 - 


i 0.2990 


0.5652 - 


2 0.2904 




2 


0.4514 - 


i 0.5074 


0.5041 - 


i 0.5071 


0.5387- 


i 0.4916 




3 


0.4141 - 


i 0.7212 


0.4715 - 


2 0.7199 


0.5051 - 


2 0.6976 



Table IV. QNFs for neutral scalar perturbations for ro = 0.1, 0.3, 0.4 for the solution of Refs. [HE]- This solution is obtained 
in the context of gravity plus nonlinear electrodynamics 4 . The solution has also been extended to the case when higher order 
curvature corrections are included in the gravitational action Ref. [S]. The parameter ro is a length scale related to the charge. 



I 
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= 0.3 


ro 


= 0.4 








0.1104- 


i 0.1150 


0.1241 - 


i 0.1104 


0.1275 


- i 0.1030 




1 


0.0954 - 


i 0.3548 


0.1073 - 


i 0.3440 


0.0993 


- i 0.3314 


1 





0.3051 - 


i 0.0991 


0.3452 - 


i 0.0989 


0.3771 


- i 0.0940 




1 


0.2779 - 


i 0.3098 


0.3227 - 


i 0.3065 


0.3549 


- i 0.2896 




2 


0.2416- 


i 0.5303 


0.2924 - 


i 0.5232 


0.3220 


- i 0.4940 




3 
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i 0.7531 


0.2550 - 


i 0.7422 


0.2806 


- i 0.7021 


2 





0.5060 - 


i 0.0980 


0.5716 - 


i 0.0982 


0.6254 


- z 0.0936 




1 


0.4871 - 


i 0.2990 


0.5564 - 


i 0.2984 
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- i 0.2834 




2 


0.4573 - 


i 0.5083 


0.5321 - 


i 0.5054 


0.5885 


- 2 0.4783 




3 


0.4207 - 


i 0.7224 


0.5023 - 


i 0.7170 


0.5583 


- i 0.6775 



Table V. QNFs for neutral scalar perturbations for ro = 0.1, 0.3, 0.4 and m — 1 for the model of Ref. [6]. 
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i 
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0.1089 


- % 0.1123 
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i 0.0970 




1 
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i 
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0.0940 
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1 
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i 


0.0980 


0.3009 
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1 
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i 


0.3074 


0.2742 


- i 0.3034 


0.3153- 


i 0.2667 




2 
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i 


0.5268 
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0.2737 - 


i 0.4581 




3 


0.1737- 


i 


0.7486 
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- i 0.7378 
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i 0.6545 


2 
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i 


0.0968 


0.4990 


- i 0.0960 


0.5689 - 


i 0.0858 




1 


0.4631 - 


i 


0.2958 
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- i 0.2929 


0.5531 - 


i 0.2599 




2 


0.4316 - 


i 


0.5034 


0.4515 


- i 0.4979 


0.5248 - 


i 0.4393 




3 


0.3925 - 


i 


0.7158 


0.4158 


- i 0.7076 


0.4871 - 


i 0.6236 



Table VI. QNFs for neutral scalar perturbations for the model of Ref. [?]■ The 5 = values reproduce those of Ref. [15] for the 
Schwarzschild case. 
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0.1032- 
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1 
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0.2587- 
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i 0.5309 
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i 0.5431 
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i 0.7724 
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0.6601 - 
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0.4584 - 


i 0.2976 


0.4451 - 


i 0.3031 




2 


0.6345 - 


i 0.4942 


0.4259 - 


i 0.5068 


0.4095 - 


z 0.5166 




3 


0.6018- 


% 0.7012 


0.3854 - 


i 0.7208 


0.3648 - 


2 0.7351 



Table VII. QNFs for neutral scalar perturbations for the solution of Ref. [Hj obtained for a system of gravity coupled to a 
phantom scalar field. The parameters have been fixed in order to normalize to unity the black hole mass, c = — 37r/(26) and 
Po — 3 andb has been set to b — 0.1, 0.5, 1. 



Appendix B: Tables of QNFs for charged scalar perturbations 
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1 


0.3005 - 


z 0.3144 


0.3369 - 


z 0.3189 


0.3729 - 


i 0.3200 




2 


0.2540 - 


z 0.5292 


0.2517- 


i 0.5192 


0.2479 - 


z 0.4841 




3 


0.1635 - 


z 0.7248 


0.0923 - 


i 0.6242 


0.0938 - 


i 0.5115 


2 





0.5182- 


z 0.0987 


0.5564- 


z 0.1002 


0.5987- 


z 0.1012 




1 


0.5000 - 


z 0.3012 


0.5397 - 


i 0.3051 


0.5831 - 


z 0.3071 




2 


0.4708 - 


z 0.5106 


0.5096 - 


z 0.5151 


0.5488 - 


z 0.5167 




3 


0.4323 - 


i 0.7219 


0.4607- 


j 0.7238 


0.4791 - 


z 0.7197 



Table VIII. QNFs for charged scalar perturbations and for the solution of Ref. [HE]. We normalized m and e to unity. 
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0.1447- 


z 0.1201 


0.1895 - 


i 0.1247 


0.2383 - 


i 0.1251 




1 


0.0724 - 


z 0.3463 


0.1684 - 


z 0.3435 


0.1368- 


i 0.2724 


1 





0.3266 - 


z 0.1013 


0.3647 - 


z 0.1044 


0.4060 - 


i 0.1072 




1 


0.3002 - 


z 0.3150 


0.3355 - 


i 0.3216 


0.3690 - 


z 0.3266 




2 


0.2535 - 


i 0.5303 


0.2500 - 


z 0.5222 


0.2407 - 


z 0.4907 




3 


0.1626- 


z 0.7262 


0.0926 - 


z 0.6279 


0.0871 - 


2 0.5162 


2 





0.5181 - 


z 0.0989 


0.5559 - 


z 0.1010 


0.5972 - 


z 0.1031 




1 


0.4998 - 


z 0.3017 


0.5388 - 


z 0.3075 


0.5807- 


z 0.3132 




2 


0.4705 - 


i 0.5116 


0.5081 - 


z 0.5193 


0.5448 - 


z 0.5268 




3 


0.4317- 


z 0.7234 


0.4581 - 


z 0.7295 


0.4721 - 


2 0.7322 



Table IX. QNFs for charged scalar perturbations and for the solution of Ref. [B]. We normalized m and e to unity. 
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0.1459- 


z 0.1199 


0.1929 


- i 0.1224 


0.2473 - 


z 0.1215 
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0.0874 - 


z 0.3515 


0.0526 


- z 0.3059 


0.0703 - 


i 0.2557 
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0.3275 - 


z 0.1012 


0.3689 


-i 0.1038 


0.4175 - 


i 0.1054 




1 


0.3014- 


z 0.3145 


0.3407 


- z 0.3194 


0.3828- 


z 0.3227 




2 


0.2559- 


z 0.5292 


0.2573 


-z 0.5199 


0.2462 - 


z 0.4930 




3 


0.1704- 


z 0.7238 


0.1008 


- z 0.6282 


0.0830 - 


i 0.5105 


2 





0.5192- 


z 0.0988 


0.5609 


- z 0.1004 


0.6103- 


i 0.1015 




1 


0.5010- 


z 0.3013 


0.5443 


- z 0.3056 


0.5951 - 


z 0.3080 




2 


0.4719 - 


z 0.5108 


0.5145 


-z 0.5160 


0.5620 - 


z 0.5184 




3 


0.4335 - 


z 0.7222 


0.4660 


- z 0.7249 


0.4943 - 


i 0.7257 



Table X. QNFs for charged scalar perturbations and for the solution of Ref. [7]. We normalized m and e to unity. 
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